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Abstract 

In this series of eight papers we present the applications of 
methods from wavelet analysis to polynomial approxima- 
tions for a number of accelerator physics problems. In this 
part we consider orbital motion in transverse plane for a 
single particle in a circular magnetic lattice in case when 
we take into account multipolar expansion up to an arbi- 
trary finite number. We reduce initial dynamical problem to 
the finite number (equal to the number of n-poles) of stan- 
dard algebraical problem and represent all dynamical vari- 
ables via an expansion in the base of periodical wavelets. 



1 INTRODUCTION 

This is the second part of our eight presentations in which 
we consider applications of methods from wavelet analy- 
sis to nonlinear accelerator physics problems. This is a 
continuation of our results from [l]-[8], which is based on 
our approach to investigation of nonlinear problems - gen- 
eral, with additional structures (Hamiltonian, symplectic or 
quasicomplex), chaotic, quasiclassical, quantum, which are 
considered in the framework of local (nonlinear) Fourier 
analysis, or wavelet analysis. Wavelet analysis is a rela- 
tively novel set of mathematical methods, which gives us a 
possibility to work with well-localized bases in functional 
spaces and with the general type of operators (differential, 
integral, pseudodifferential) in such bases. In this part we 
consider orbital motion in transverse plane for a single par- 
ticle in a circular magnetic lattice in case when we take into 
account multipolar expansion up to an arbitrary finite num- 
ber. We reduce initial dynamical problem to the finite num- 
ber (equal to the number of n-poles) of standard algebraical 
problem and represent all dynamical variables as expansion 
in the base of periodical wavelet functions. Our consider- 
ation is based on generalization of variational wavelet ap- 
proach from part 1 . After introducing our starting points 
related to multiresolution in section 3, we consider meth- 
ods which allow us to construct wavelet representation for 
solution in periodic case in section 4. 
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2 PARTICLE IN THE MULTIPOLAR 
FIELD 

The magnetic vector potential of a magnet with 2n poles in 
Cartesian coordinates is 



A = ^K n f n {x,y), 



(1) 



where /„ is a homogeneous function of x and y of order n. 
The real and imaginary parts of binomial expansion of 



f n (x,y) = [x + iyY 



(2) 



correspond to regular and skew multipoles. The cases 
n = 2 to n = 5 correspond to low-order multipoles: 
quadrupole, sextupole, octupole, decapole. Then we have 
in particular case the following equations of motion for sin- 
gle particle in a circular magnetic lattice in the transverse 
plane (x, y) ([9] for designation): 
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and the corresponding Hamiltonian: 
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Then we may take into account arbitrary but finite num- 
ber in expansion of RHS of Hamiltonian (Q) and from our 
point of view the corresponding Hamiltonian equations of 
motions are not more than nonlinear ordinary differential 
equations with polynomial nonlinearities and variable co- 
efficients. 



3 WAVELET FRAMEWORK 

Our constructions are based on multiresolution approach. 
Because affine group of translation and dilations is inside 
the approach, this method resembles the action of a mi- 
croscope. We have contribution to final result from each 
scale of resolution from the whole infinite scale of spaces. 
More exactly, the closed subspace Vj(j G Z) corresponds 
to level j of resolution, or to scale j. We consider a r-regular 
multiresolution analysis (MRA) of L 2 (R") (of course, we 
may consider any different functional space) which is a se- 
quence of increasing closed subspaces Vf 

...V-2 c V_i cVoCVicVjC ... (5) 
satisfying the following properties: 

f]Vj = 0, ~\JVj = L 2 (R"), 
jez jez 

f(x)eVj <=> f(2x)eV j+1 , 
f(x) e V <=> f(x -k)e V Q , Vfc g Z". (6) 

There exists a function ip G V) such that {(p ^{x) = (p(x — 
k), k G Z"} forms a Riesz basis for Vq. The function 
(p is regular and localized: tp is C r_1 , <^ r_1 ) is almost 
everywhere differentiable and for almost every x G R™, 
for every integer a < r and for all integer p there exists 
constant C p such that 

\d a <p(x)\<C p (l + \x\)- p (7) 

Let <p(x) be a scaling function, ip(x) is a wavelet function 
and ifi(x) — ip(x — i). Scaling relations that define <p, if) 
are 
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Let indices £, j represent translation and scaling, respec- 
tively and 

tpjl(x) = 2 j / 2 ip(2 j x-£) (10) 

then the set {<fij,k}, k G Z™ forms a Riesz basis for Vj. The 
wavelet function tp is used to encode the details between 
two successive levels of approximation. Let Wj be the or- 
thonormal complement of Vj with respect to Vj+i: 

V 3+ i=V 3 ^W r (11) 

Then just as Vj is spanned by dilation and translations of 
the scaling function, so are Wj spanned by translations and 
dilation of the mother wavelet ipjk (x), where 

ipjk{x) = 2 j / 2 ip(2 j x- k). (12) 



All expansions which we used are based on the following 
properties: 

{i>jk}, j, k G Z is a Hilbertian basis of £ 2 (R) 
Wjk}j>o.kez is an orthonormal basis forL 2 (R), 

oo 

L 2 (H) = V ^W J , (13) 

j=o 

or {ifio,k,ipj,k}j>o,kez 

is an orthonormal basis fori 2 (R) . 

Fig.l and Fig. 2 give the representation of some function 
and corresponding MRA on each level of resolution. 




Figure 2: MRA representation. 



4 VARIATIONAL WAVELET APPROACH 
FOR PERIODIC TRAJECTORIES 

We start with extension of our approach from part 1 to 
the case of periodic trajectories. The equations of motion 
corresponding to Hamiltonian (4) may also be formulated 
as a particular case of the general system of ordinary dif- 
ferential equations dxi/dt = fi(xj,t), = 1, ...,n), 
< t < 1, where /j are not more than polynomial func- 



tions of dynamical variables Xj and have arbitrary depen- 
dence of time but with periodic boundary conditions. Ac- 
cording to our variational approach from part 1 we have the 
solution in the following form 

x i {t)=x i {0)+J2*i<Pk(t), x i (0)=x i {l), (14) 

where are again the roots of reduced algebraical sys- 
tems of equations with the same degree of nonlinearity and 
ipk (t) corresponds to useful type of wavelet bases (frames). 
It should be noted that coefficients of reduced algebraical 
system are the solutions of additional linear problem and 
also depend on particular type of wavelet construction and 
type of bases. This linear problem is our second reduced 
algebraical problem. We need to find in general situation 
objects 
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but now in the case of periodic boundary conditions. Now 
we consider the procedure of their calculations in case 
of periodic boundary conditions in the base of periodic 
wavelet functions on the interval [0,1] and corresponding 
expansion (14) inside our variational approach. Periodiza- 
tion procedure gives us 

0j,k( X ) = ^2f],k(x ~ i) (16) 

eez 

4>j,k( x ) = y,v>i,fc(ff - 1) 

So, <p,ip are periodic functions on the interval [0,1]. Be- 
cause (fj_k = fj,k' if k — fc'mod^), we may consider 
only < k < V and as consequence our multiresolution 
has the form (J % = L 2 [0, 1] with % = span^fc}^" 1 

[10]. Integration by parts and periodicity gives useful re- 
lations between objects (15) in particular quadratic case 
(d = di+d 2 ): 



A 



di ,d 2 
fci ,ko 



A 



0,(2 

0,fc 2 -fc, 



\ d 



So, any 2-tuple can be represent by A^. Then our sec- 
ond additional linear problem is reduced to the eigenvalue 
problem for {A^}g <A .< 2 i by creating a system of 2- 5 ho- 
mogeneous relations in A d and inhomogeneous equations. 
So, if we have dilation equation in the form <p(x) = 
\/2 ^2 keZ hkf(2x — fc), then we have the following ho- 
mogeneous relations 
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or in such form AX d = 2 d X d , where X d = {Af} < fe < 2 i. 
Inhomogeneous equations are: 



J2 MfAf = d\2 



-Hi 



(18) 



where objects M^(|^| < N — 2) can be computed by re- 
cursive procedure 



(19) 



M\ =< x k ,^ £ >= E ( ) n k ~ j Mi, M* = 1. 

3=0 

So, we reduced our last problem to standard linear alge- 
braical problem. Then we use the same methods as in part 
1. As a result we obtained for closed trajectories of orbital 
dynamics described by Hamiltonian (4) the explicit time 
solution (14) in the base of periodized wavelets (16). 
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